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Abstract 

Missions to rendezvous with or capture an asteroid present signifi¬ 
cant interest both from a geophysical and safety point of view. They 
are key to the understanding of our solar system are as well step¬ 
ping stones for interplanetary human flight. In this paper, we focus 
on a rendezvous mission with 2OO6RH120 , an asteroid classified as a 
Temporarily Captured Orbiter (TCO). TCOs form a new population 
of near Earth objects presenting many advantages toward that goal. 

Prior to the mission, we consider the spacecraft hibernating on a Halo 
orbit around the Earth-Moon’s libration point. The objective is 
to design a transfer for the spacecraft from the parking orbit to ren¬ 
dezvous with 2OO6RH120 while minimizing the fuel consumption. Our 
transfers use indirect methods, based on the Pontryagin Maximum 
Principle, combined with continuation techniques and a direct method 
to address the sensitivity of the initialization. We demonstrate that 
a rendezvous mission with 2OO6RH120 can be accomplished with low 
delta-v. This exploratory work can be seen as a first step to identify 
good candidates for a rendezvous on a given TCO trajectory. 

Keywords: Temporarily Captured Objects, Three-body Problem, Opti¬ 
mal control, Indirect Numerical Methods 
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1 Introduction 


The motivation for our work is to study asteroid capture missions for a spe¬ 
cific population of near Earth objects. The targets, Temporarily Captured 
Orbiters (TCO), are small asteroids that become temporarily-captured on 
geocentric orbits in the Earth-Moon system. They are characterized as sat¬ 
isfying the following constraints: 

• the planetocentric Keplerian energy Epianet < 0; 

• the planetocentric distance is less than three Earth’s Hill radii (e.g., 
3Rh,(b ~ 0.03 AU); 

• it makes at least one full revolution around the Earth in the Earth-Sun 
co-rotating frame, while satisfying the first two constraints. 

In regard to the design of a round trip mission, the main advantage of the 
TCOs lies in the fact that those objects have been naturally redirected to 
orbit the Earth. This contrasts with recently proposed scenarios to design, 
for instance, a robotic capture mission for a small near-Earth asteroid and 
redirect it to a stable orbit in the EM-system, to allow for astronaut visits 
and exploration (e.g. the Asteroid Redirect Mission (ARM)). 

In this paper we focus on 2 OO 6 RH 120 , a few meters diameter near Earth 
asteroid, officialy classified as a TCO. 2 OO 6 RH 120 was discovered by the 
Catalina Sky Survey on September 2006. Its orbit from June 1st 2006 to 
July 31st 2007 can be seen on Figure generated using the Jet Propulsion 
Laboratory’s HORIZONS database which gives ephemerides for solar-system 
bodies. The period June 2006 to July 2007 was chosen to include the portion 
of the orbit within the Earth’s Hill sphere with a margin of about one month. 
We can also observe that 2 OO 6 RH 120 comes as close as 0.72 normalized units 
from Earth-Moon barycenter. In [H], the authors investigate a population 
statistic for TCOs. Their work is centered on the integration of the tra¬ 
jectories for 10 million test-particles in space, in order to classify which of 
those become temporarily-captured by the Earth’s gravitational field - over 
eighteen-thousand of which do so. This suggests that 2 OO 6 RH 120 is not the 
only TCO and that it is relevant to compute a rendezvous mission to this 
specific asteroid to get insight whether TCOs can be regarded as possible 
targets for transfers with small fuel consumption, and thus cost. 

In this paper, our goal is to solve the optimal control problem that consists 
in performing a rendezvous mission from a Halo orbit around the Earth- 
Moon libration point L 2 to the 2 OO 6 RH 120 asteroid, while minimizing the fuel 
consumption. The choice of targets for our rendezvous mission sets us apart 
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Figure 1; Orbit of 2 OO 6 RH 120 in the Earth-Moon CR3BP rotating (left) and 
inertial (right) reference frame. 


from the existing literature where transfers are typically designed between 
elliptic orbits in the Earth-Moon or other systems |T9] and | 20 ], or to a 
Libration point HB m\ and [2H 1^ . Rendezvous missions to asteroids in 
the inner solar system can be found in |101 [T7| but they concern asteroids on 
elliptic orbits which is not the case for us since TCOs are presenting complex 
orbits and therefore require a different methodology. Our assumption on 
the hibernating location for the spacecraft, a Halo orbit around the Earth- 
Moon unstable Libration points L 2 , is motivated in part from the successful 
Artemis mission [221123] and in part from the constraint on the duration of 
the mission, mostly impacted by the time of detection of the asteroid. Indeed, 
the Arthemis mission demonstrated low delta-v (An) station keeping on Halo 
orbits around Li and L 2 . While this first study focuses on similar Halo 
orbits to the Artemis mission, the study needs to be extended to analyze a 
variation of Halo orbits both in shape (to include the eight-shaped Lissajous 
orbits) and also in 2 ;-excursion or energy. As a hrst approach, in this work 
we search the best departure point on the prescribed Halo orbit without 
practical consideration related to synchronization with the asteroid orbit. 
The primary objective for the mission is to maximize the hnal mass, but 
since the Earth capture duration of the asteroid is limited it imposes a time 
constraint which need to be addressed eventually. In this paper however, we 
consider rendezvous missions possibly over long period of time to determine 
what is the correlation between the required Av and the duration of the 
mission. A forthcoming paper will bridge our results to the TCOs detection 
aspect. 

Many methodologies have been developed over the past decades to design 
optimal transfers in various scenarios. Due to the complexity of the TCO 


3 










orbits and the nature of the mission, techniques based on analytical solutions 
such as in [16] for circular Earth orbits are not suitable and we use a numerical 
approach. A survey on numerical methods can be found in |H], and for 
reasons related to the specihcs of our problem we choose to use a deterministic 
approach based on tools from geometric optimal control versus an heuristic 
method such as in mu and |2Zl l29] . More precisely, we use an indirect 
method based on the maximum principle, as well as a direct method and 
continuation techniques to address the difficulties of initialization for our 
numerical scheme. Additionally, we fix the structure of the control norm to be 
a piecewise constant function with the magnitude of the thrust either zero or 
maximum with at most three maximum magnitude arcs. Our techniques are 
illustrated on the 2 OO 6 RH 120 TOO and four others coming from the database 
of [13]. Validation of our approach can be seen by comparing the work to 
ra. in which the authors develop a low Av asteroid rendezvous mission 
that makes use of a Halo orbit around Earth-Moon L 2 - Their situation is, 
however, different from ours, in that they have carefully chosen a specihc 
asteroid for rendezvous. With a one-year transfer time, the Av value they 
realize is 432 m/s, which is comparable to the Av values presented here in a 
less-ideal scenario. 


2 Optimal Control Problem and Numerical Al¬ 
gorithm 

2.1 Model 

In this paper, the circular restricted three-body problem |18j is used to ap¬ 
proximate the spacecraft dynamics. This is justihed by the fact that a TCO 
can be assumed of negligible mass, and that the spacecraft evolving in the 
TCO’s temporary capture space is therefore attracted mainly by two primary 
bodies, the Earth and the Moon. 

The CR3BP model is well known, let us however go through some nota¬ 
tions fo the sequel of the paper. We denote by {x(t),y{t), z(t)) the spatial 
position of the spacecraft at time t. In the rotating coordinates system, and 
under proper normalization, the primary planet identihed here to the Earth, 
has mass mi = 1 — y and is located at the point {—fi, 0,0); while the second 
primary, identihed to the Moon, has a mass of m 2 = /i and is located at 
(1 — /i, 0,0), where fi = 1.2153e — 2. The distances of the spacecraft with 
respect to the two primary bodies are given by pi = -s/{x yY + 2/^ + 2 ;^, 
Q 2 = \/(x — 1 -f yY + y'^ + respectively. The potential and kinetic energy, 
respectively V and K of the system are given by 
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x'^ + y‘^ 1-H /i /i(l - /i) 

2 ^2 2 ’ 


• 2 I • 2 I • 2\ 

-(a: +y + z ). 


( 1 ) 


We assume a propulsion system for the spacecraft is modeled by adding terms 
to the equations of motion depending on the thrust magnitude and some 
parameters related to the spacecraft design. The mass of the spacecraft is 
denoted by m. Under those assumptions, we have the following equations of 
motion: 


.. dV ^ T^ax 

x-2y = —^ -Ml, 

ox m 


oy m 


z = 


dV Taiax 
-T- +-Ms (2) 


dz 


m 


where u{.) = n 2 (-),M 3 (-)) is the control, and satisfies the constraint 

IImII = ^Ju\ + u^ + u\ < 1, and with Tma^ /m normalized. A first integral of 
the free motion is given by the energy of the system E = K — V. We will 
later use this energy value to analyze the choice of the rendezvous point and 
the parking orbit for the spacecraft. It is well know that the uncontrolled 
motion of the dynamical system has five equilibrium points defined as the 
critical points of the potential V. Three of them Li,L 2 and L 3 are aligned 
with the Earth-Moon axis and have been shown to be unstable, while the 
two others are stable and are positioned to form equilateral triangles in the 
plane of orbit with the two primaries. Since our goal is to maximize the final 
mass we must include the differential equation governing the variation of the 
mass along the transfer: 

m = -/3Tmax||M|| (3) 

where the parameter f3, the thruster characteristic of our spacecraft, is given 
hj (5 = (it is the inverse of the ejection velocity Ue), with Ap the specific 
impulse of the thruster and go the acceleration of gravity at Earth sea level. 


2.2 Optimal Control Problem 


Let q = {qs,qvY where = {x,y,zy represents the position variables and 
q^ = {x,y,zy the velocity ones. From section 2.1, 
an affine control system of the form: 


our dynamical system is 


q = Fo{q) + Fi{q)ui (4) 

i=\ 

where the drift is given by Fo(g) = (x, y, z, 2y+x— _ mU-i+z^) ^ _2xy 

B\ 82 

y ~ ~ ~ ~ control vector fields are Tj(g) = e^+i 

with Cj forming the orthonormal basis of IR®. 
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We consider the rendezvous transfer from an initial point on 

a parking orbit Oq G IR® to a final position and velocity on the 

2 OO 6 RH 120 orbit. Note that the initial and hnal positions and velocities are 
variables of the global optimization problem. The criterion to maximize is 
the hnal mass which is equivalent to minimizing the fuel consumption or the 
An = ^i'^ce the mass evolves proportionally to the norm of 

the thrust, our criterion is equivalent to the minimization of the Li-norm of 
the control; 

min [ \\u{t)\\dt, (5) 

“ 6 " JtQ 


where U = {«(.);measurable bounded and ||M(t)|| < 1 for almost all t}, Iq 
and tf are respectively the initial and hnal time. Remark that since we choose 
Oo to be a Halo orbit around a libration point, it is uniquely determined by 
a single point of the orbit using the uncontrolled CR3BP dynamics. This 
fact will play an important role for one of the necessary optimality condi¬ 
tions below. The large number of variables in our problem signihcantly adds 
complexity to the search for a solution. In particular, in the case of free hnal 
time we expect an inhnite time horizon with a control structure that mimics 
impulsive maneuvers. To simplify our optimal problem we have two options, 
either we hx the transfer time or we hx the structure of the control. If we 
hx the hnal time, the sensitivity of the shooting method can be addressed by 
using the solution of a smoother criterion than the Li-norm, for instance the 
L 2 -norm, and linking it to the target criterion by a continuation procedure, 
see na. In this paper, however, we take a diherent approach and decide to 
hx the structure of the control. In the sequel we focus on designing transfers 
associated to controls with a piecewise constant norm with value in { 0 , 1 } 
and four switchings. More precisely, we consider control functions that are 
piecewise continuous such that there exists times to < ti < t 2 < < tf 

with 


f 1 if t e {to]ti)U {t 2 ;t 3 )U {tiitf) 
\ 0 if G (ti;^ 2 ) U (^3!^ 4 )- 


( 6 ) 


Here the hnal time t/ is free. Note that our numerical method will be able 
to select the best control strategy even if it has less than three boosts. 


2.3 Necessary Conditions for Optimality 

The maximum principle provides hrst order necessary conditions for a tra¬ 
jectory to be optimal j^B]- Details regarding the application of the maxi¬ 
mum principle to orbital transfers can be found in many references, includ¬ 
ing jH [26]. We denote by X{t) = {q{t),m{t)) G IR®^^ the state, where 
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Figure 2: Control function modeling thrust impulses over time. 


q = {x, y, z, X, y, z) is the position and velocity of the spacecraft with m 
its mass. The maximum principle applied to our optimal control prob¬ 
lem states that if (q'(-),m(-),«(■)) = (X(-),m(-)) is an optimal solution de- 
hned on then there exists an absolutely continuous adjoint state 

(p°,Px(-)) = dehned on [to,tf] such that: 

• {p^yPxi')) 7 ^ 0, Vf G and < 0 is a constant. 

• Let H, the Hamiltonian, be H{t,X{t),p^,px{t),u{t)) = p°||u(f)|| -|- 
{pxit),X{t)), then 

dH 

x{t) = ■^{t,X{t),p°,px{t),u{t)), for a.e. t e [to,tf], (7) 

dH 

Px{t) = --^{t,X{t),p^,px{t),u{t)), for a.e. t e [to,tf], (8) 
where (,) denotes the inner product. 

• H{t,X{t),p°,px{t),u{t)) = max\\^\\<iH{t,X{t),p°,px{t),iy), Vfs.t. ||M(t)|| 
1 (maximization condition). 

• '1/(L) = 0 for i = 1, • • • , 4. 

• H{tf,X{tf),p°,pxitf),u{tf)) = 0, if tf is free. 

• {Pqito)y = 0 (initial transversality condition). 

• Pmitf) = 0. 

The function T(-) is the so-called switching function corresponding to the 
problem with an unrestricted control strategy and we have 4/(f) = Po + 

rp ( Ilpdbll _ ^ (t)B\ 
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The maximization condition of the Hamiltonian H is used to compute 
the control on U [t 2 ;^ 3 ] U [^ 4 ;^/] and we have u{t) = for all t G 

[fo;ti] U [f 2 ;^ 3 ] U [ti;tf]. The initial transversality condition reflects the fact 
that the initial departing point is free on the Halo orbit Oq. Remark that 
since the data for the TCO’s trajectory are given as ephemerides, there are no 
dynamics equations to describe those orbits in the CR3BP. Thus, we cannot 
compute the tangent space to a TCO point and we cannot extract a transver¬ 
sality condition for Pq{tf) at the rendezvous point. Since we expect numerous 
local extrema for this optimal control problem, it is however preferable to 
solve the problem for fixed rendezvous points on a discretization of the TCO 
orbit. 

2.4 Shooting Method 

Our numerical method is based on the necessary conditions given in section 


expressed using the maximization condition. Then, we have Z{t) = <h(Z(f)) 
where $ comes from equations Q and ([^. The goal is to find Z(to), ti, ^ 2 , 
^ 3 , ^4 and tf such that the following conditions are fulfilled; 

1. = 0 for z = 1, • • ■ 4; 

2. X{tf) is the prescribed rendezvous point; 

3. X(fo) G Oq and the initial transversality condition is verified; 

4. H(tf) = 0. 

The problem has been transformed into solving a multiple points boundary 
value problem. More specifically, we must find the solution of a nonlinear 
equation S'(Z(to),h,^ 2 ,^ 3 ,^ 4 ,^/) = 0 , where S is usually called the shooting 
function. The evaluation of the shooting function is performed using the 
high order numerical integrator DOP853, see na. The search for a zero of 
the shooting function is done with the quasi-Newton solver HYBRD of the 
Fortran minpack package. Since S{.) is nonlinear, the Newton-like method 
is very sensitive to the initial guess. This leads us to consider heuristic ini¬ 
tialization procedures. We combine two types of techniques, a direct method 
and a continuation method. The discretization of the TCO’s orbit requires 
the study of thousands of transfers, we use a direct method for a dozen of 
rendezvous and expand to other points on the orbit using a continuation 
scheme. The motivation is that direct methods are very robust but time 
consuming while continuation method succeeded for our problem in most 


2.3 Let Z{t) = {X{t),px{t)), t G [to]tf], and u{q,p) the feedback control 
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cases very rapidly. The direct method uses the modeling language Ampl, 
see [l2], and the optimization solver IpOpt, see j28], with a second order ex¬ 
plicit Runge-Kutta scheme. Details on advanced continuation methods can 
be found in [13] and |5], but we use a discrete continuation which is enough 
for our needs. 

Moreover, note that in order to fulfill the initial transversality condition, 
we first prescribe X(to) on the parking orbit and find a zero of the shooting 
function satisfying all the other conditions. Afterward, we do a continuation 
on X(fo) along the departing parking orbit in the direction that increases the 
hnal mass. Once a framing of the best X{to) with respect to the final mass is 
found, we perform a final single shooting to satisfy the initial transversality 
condition along with the other conditions. This decoupling in the search of an 
extremal is motivated by the fact that we could very likely find a local maxi¬ 
mum on X(fo) rather than a local minimum because of the periodicity of the 
initial parking orbit. We avoid this fact by first manually ensuring that the 
X (to) we find will be the one for the best final mass and not the worst. How¬ 
ever, our continuation procedure on X(to) does not always succeed, mainly 
because of the high nonlinearity of the shooting function, as the trajectories 
we find can be very long. Even if some of the extremals we find do not satisfy 
the initial transversality condition with the aimed accuracy (typically a zero 
of the shooting function is deemed acceptable if ||5'(Z(fo), 2,3,4, /)|| < 10 “®), 

they are still rather close to satisfy it (of the order of 10 “^). 

3 Results 

3.1 Rendezvous to 2 OO 6 RH 120 

The objective of this section is to present an analysis of the evolution of the 
fuel consumption with respect to the location of the rendezvous point on the 
2 OO 6 RH 120 orbit. As mentioned before, we restrict the study to transfers 
with at most three boosts. The primary goal is to obtain insights on the 
variation of fuel consumption based on the features of the rendezvous point 
such as its distance from L 2 , or its energy. The spacrecraft characteristics are 
assumed to be an initial mass of 350 kg, a specific impulse Ap of 230 s. and 
a maximum thrust T^ax of 22 N. The Halo orbit from which the spacecraft 
is departing is chosen to have a 2 ;—excursion of 5000 km around the EM 
libration point L 2 , see Figure]^ The choice for this specific Halo orbit is mo¬ 
tivated by the successful Artemis mission [221123] that used Halo orbits with 
similar characteristics. The point corresponding to the positive z—excursion 
is (?HaioL 2 = (1.119,0,0.013, 0., 0.180,0), and the period of this particular 
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Halo orbit is tHaioL 2 = 3.413 in normalized time units or approximatively 
14.84 days. 



y(LD) 


Figure 3; Halo orbit from which the spacecraft is departing, z—excursion of 
5000 km around the EM libration point L 2 . 

During the period represented in Figure [T| asteroid 2 OO 6 RH 120 does 17 
clockwise revolutions around the origin of the CR3BP frame, and 3.6 rev¬ 
olutions in Earth inertial reference frame. The evolution of the energy of 
2 OO 6 RH 120 and its distance to the Earth-Moon libration point L 2 as the as¬ 
teroid evolves on its orbit are given in Figure]^ 




Figure 4: Evolution of the energy (left) and of the Euclidean distance to the 
L 2 Libration point (right) for 2 OO 6 RH 120 • For energy, the horizontal line 
represents the energy of the Halo periodic orbit around L 2 which is about 
-1.58. 

To analyze the variations of the fuel consumption with respect to the 
rendezvous point on the orbit, we discretize uniformly the orbit of 2 OO 6 RH 120 
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using 6 hours steps. For each rendezvous point of this discretization, we 
compute an extremal transfer (i.e. a solution of the maximum principle) 
with free final time using the techniques explained in section Figure 
shows the evolution of the final mass, the An and of the duration of the 
transfer with respect to the rendezvous point on the 2 OO 6 RH 120 orbit for a 
spacecraft departing from the Halo L 2 orbit and corresponding to the three 
boosts control strategy. As explained in section 2 H, some of the departure 
point on the Halo orbit are not fully optimized - this is the case for about 
two thirds of them. Also note that a departure point different from gHaioL 2 
implies a drift phase whose duration is not included in the transfer duration. 
It can be observed that the final mass has many local extrema and that 
the variation of the duration of the mission is not continuous (contrary to 
what we would expect). It is most likely due to local minima or to the fact 
that the value function of the problem can present discontinuities. Figure 





Figure 5: Evolution of the final mass (top left), of the An (top right) and 
transfer duration (bottom) from the Halo orbit around L 2 with respect to 
rendezvous points on the orbit of 2 OO 6 RH 120 • 

shows the evolution of the departure point on the hibernating orbit for the 
spacecraft with respect to the rendezvous point on 2 OO 6 RH 120 (left) as well 
as the three most frequent departure points on the initial Halo orbit (right). 
More precisely, we represents the optimized argument of (up to the 
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quadrant: arctan{z{^)/y{^))). Note that the initial position on the Halo 
orbit directly depends on the optimized initial drift time. Since this initial 
drift phase has not always been successfully optimized, this figure has to be 
interpreted with caution. However, we can see that the departure position 
on the initial periodic orbit seems to always be close to a multiple of vr. 




tRH120 (days) 


Figure 6 : The left picture represents the evolution of (?/, ^)(0) argument with 
respect to time of rendezvous on 2 OO 6 RH 120 • The right picture represents 
the three most frequent departure points on the initial Halo orbit. 

Comparing the evolution of the final mass from Figure and the evolution 
of 2 OO 6 RH 120 energy from Figure we can see that the best final masses 
are obtained on the first half of 2 OO 6 RH 120 trajectory, that is for rendezvous 
points which energies are closer to the departing energy. For the best transfer, 
the energy difference between the rendezvous point on 2 OO 6 RH 120 and the 
departing point on the Halo orbit is about 0.046. The rendezvous point 
with the closest energy to the initial orbit is only slightly before the optimal 
rendezvous point. Its final mass is 319.3 kg which is only 0.5 kg worst than 
the best final mass. This remark suggests strongly that a small difference in 
energy between the rendezvous point on 2 OO 6 RH 120 orbit and the departing 
point for the spacecraft on the Halo orbit is advantageous. 

The best transfer is represented on Figures and some relevant data 
is also presented in Table This transfer has a An of 203.6 m/s. The 
rendezvous takes place on June 26*^ 2006 and lasts 415.5 days which would 
require to detect and launch the mission about 14 months before June R* 
2006. On Figure we display the orbit of 2 OO 6 RH 120 in both the rotating 
and inertial frame with the rendezvous point for the best transfer. The best 
transfer exhibits 14 revolutions around the origin (in the rotating frame) and 
has a significant variation in the z—coordinate with respect to the EM plane. 
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Table [T| Best Transfer Data from Halo-L 2 


Parameter 

Symbol 

Value 

Transfer Duration 

tf 

415.8 days 

Final Mass 

mj 

319.8 kg 

Delta-v 

Av 

203.6 m/s 

Final Position 


(2.25, 3.21, -1.04) 

Final Velocities 


(2.92, -2.02, 0.46) 

Max Distance from L 2 

^max 

^L2 

10.63 LD 


Table 1: Data for the best transfer from gHaioL 2 to asteroid 2 OO 6 RH 120 • 



- Beghnlnacrflfaratar 



Figure 7: Best 3 boosts rendezvous transfer to 2 OO 6 RH 120 from a Halo orbit 
around L 2 . Inertial Frame (left). Rotating Frame (Right). Bottom: zoom 
on the start of the departure from the Halo orbit (rotating frame). 


In particular, the coordinate of the rendezvous point is —1.04 normalized 
units, that is about 400 thousand km, and the maximum coordinate along 
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Figure 8 : Orbit of 2 OO 6 RH 120 in rotating (right) and inertial (left) reference 
frame with the rendezvous point corresponding to the best transfer. 


the trajectory is 5.34 normalized units, that is about 2 million km. The 
departure point on the parking orbit occurs 4.5 days after gHaioL 2 - The control 
strategy consists of a hrst boost of 19.7 s, a second boost starting after 
154.7 days and lasting 51.1 min , hnally, the last boost starting 261.1 days 
after the second boost and lasts 13.8 s, see Figure Despite a seemingly 
short initial boost, this trajectory does not exhibit the prohle of an initial 
jump on an unstable invariant manifold. However, this trajectory exploits 
the fact that a small initial boost leads to a far location where the gravity 
held of the two primaries is small and where the second boost can efficiently 
aim at the rendezvous point. In particular, we expect the existence of other 
local minima with a larger hnal time, going further away from the initial and 
hnal positions and providing an even better hnal mass. Also note that this 
kind of strategy could not have been obtained if we had restricted the control 
structure to have two boosts rather than three. 

To contrast with the best transfer we represent the worst transfer in Fig¬ 
ures 10 This transfer has a An of 962.5 m/s (hnal mass is m/ = 22%.‘ikg) and 


a transfer time of 37.9 days. The hnal position is = (—0.65,0.35,0.92) 
and the hnal velocity is = (—0.61,1.21, —0.28). The energy at the worst 
rendezvous point is about —0.168, and the maximum distance to L 2 is 6.21 
LD. 


4 Conclusion and Future Work 

The numerical approach using an indirect method produces a low delta-v 
transfer to asteroid 2 OO 6 RH 120 departing from a Halo orbit around EM-L 2 - 
Further calculations on four synthetic TCOs obtained from the databased 
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Figure 9: The bottom graph represents the control strategy for the best 
transfer to 2 OO 6 RH 120 , and the top graph shows the evolution of the energy 
for this transfer. The energy of initial Halo orbit EHaioL 2 = -1.581360, the 
energy after 1st boost = -1.581335, which means that the first boost raised 
the energy for about 2.46e — 05. The energy after 2nd boost = -1.535022 
which is an increase of about 0.046313 and the energy at the rendezvous 
point = -1.535804 which is a decrease of about = 7.814e-04. The total energy 
difference is about 0.0456. 

produced in na suggest that similar results can be expected on a large 
sample of TCOs. Indeed, in Table we display data regarding the best 
three boosts transfer for four other TCOs. These transfers have a duration 
of about one year each and produce delta-v values between 223.9 m/s and 
344.2 m/s. The best transfer is for TCOi and occurs for an energy difference 
between terminal conhgurations of about 0.13, for TCOie the difference is 
about 0.6, 1.5 for TCOig and 2.3 for TCOn. This reinforce the idea of a 
relation between the hnal mass of the transfer and the energy difference of 
the rendezvous point with respect to the one from the Halo orbit. This will 
be explored on a larger population of TCOs with the goal to characterize 
which asteroids are better suited for a rendezvous mission with a low delta- 
V. It can also be noted that the maximum distance of the spacecraft from 
L 2 during the transfer is similar for all four transfers which indicates that 
the long drift is used to pull away from the two primaries attraction helds to 
make the second boost more efficient. 

Forthcoming work will also add the Sun as a perturbation into the model, 
such as in [H [21], to determine its impact on the transfers. Moreover, for 
practical reason it would be very interesting to consider the spacecraft parked 
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Figure 10; Worst 3 boosts rendezvous transfer to 2 OO 6 RH 120 from a Halo 
orbit around L 2 . Inertial Frame (left). Rotating Frame (Right). Bottom; 


zoom on the start of the departure from the Halo orbit (rotating frame) 

Table Best Transfer Data from Halo-L 2 for Selected TCOs 


Parameter 

Symbol 

TCOi 

TCOn 

TCO 16 

TCO 19 

Transfer Duration (days) 

tf 

362.0 

386.6 

362.2 

364.9 

Final Mass (kg) 

ruf 

316.9 

300.5 

311.0 

307.1 

Delta-v (m/s) 

Av 

223.9 

344.2 

266.1 

294.6 

Max Distance from L 2 

jmax 

^L2 

12.7 

11.5 

11.5 

12.8 

Time to (days) 

tid^n 

232.9 

196.6 

197.3 

232.3 


Table 2; Best transfer data from q'HaioL 2 to selected TCOs. 


on a Halo orbit around the ES Li libration point. 

5 Acknowledgements 

We would like to specially thank Robert Jedicke and Mikael Granvik for their 
generous access to the database of synthetic TCOs, and in general, their help 
and support for this research. Geoff Patterson and Monique Chyba were 


16 













partially supported by the National Science Foundation (NSF) Division of 
Mathematical Sciences, award ^1109937. 


References 

[1] E. Belbruno and J. Miller, Sun perturbed earth-to-moon transfers with 
ballistic capture, Journal of Guidance and Control, vol. 16, no. 4, pp. 
770-775, 1993. 

[2] C.R. Besette, D.B. Spencer, Optimal Space Trajectory Design: A 
Heuristic-Based Approach, Adavances in the Astronautical Sciences, Vol. 
124, pp.1611-1628, 2006. 

[3] J.T. Betts, Practical Methods for Optimal Control Using Nonlinear Pro¬ 
gramming. Society for Industrial and Applied Mathematics, 2001. 

[4] J.B. Caillau, B. Daoud, J. Gergaud, Minimum fuel control of the pla¬ 
nar circular restricted three-body problem. Celestial Mech. Dynam. As- 
tronom. 114, no. 1, pp. 137-150, 2012. 

[5] J.B. Caillau, O. Cots and J. Gergaud, Differential continuation for reg¬ 
ular optimal control problems. Optim. Methods Softw., 27, pp. 177-196, 
2012 . 

[6] M. Chyba, G. Patterson, G. Picot, M. Granvik, R. Jedicke, J. Vaubail- 
lon. Designing Rendezvous Missions with Mini-Moons using Geometric 
Optimal Control. Journal of Industrial and Management Optimization, 
Vol. 10, Issue 2, pp. 477-501, 2014. 

[7] M. Chyba, G. Patterson, G. Picot, M. Granvik, R. Jedicke, J. Vaubail- 
lon. Time-minimal orbital transfers to temporarily-captured natural 
Earth satellites. PROMS Series: Advances of Optimization and Con¬ 
trol With Applications, Springer Verlag, 2014. 

[8] B.A. Conway, A Survey of Methods Available for the Numerical Opti¬ 
mization of Continuous Dynamic Systems, J Optim Theory Appl 152, 
pp. 271-306, 2012. 

[9] M. Pontani, B.A. Conway, Optimal Finite-Thrust Rendezvous Trajec¬ 
tories Found via Particle Swarm Algorithm, Journal of Spacecraft and 
Rockets, Vol. 50, No. 6, November-December 2013. 


17 



[10] D.W. Dunham, R.W. Farquhar, N. Elsmont, E. Chumachenko, S. Ak¬ 
senov, Y. Fedorenko, A. Genova, J. Horsewood, R. Furfaro, J. Kidd Jr, 
Using lunar swingbys and libration-point orbits to extend human explo¬ 
ration to interplanetary destinations, 64th International Astronautical 
Congress, Beijing, China, 2013. 

[11] D.C. Folta, T.A. Pavlak, A.F. Haapala, K.C. Howell, M.A. Woodard, 
Earth-Moon libration point orbit station keeping: Theory, modeling, and 
operations, Acta Astronautica, 2013. 

[12] R. Fourer, D.M. Gay, B.W. Kernighan, AMPL: A Modeling Language 
for Mathematieal Programming, Duxbury Press, Brooks-Cole Publishing 
Company, 1993. 

[13] J. Gergaud, T. Haberkorn, P. Martinon, Low thrust minimum-fuel or¬ 
bital transfer: a homotopic appoach. Journal of Guidance, Control and 
Dynamics, Vol.27, No.6, pp. 1046-1060, 2004. 

[14] M. Granvik, J. Vauballion, R. Jedicke, The population of natural Earth 
satellites, Icarus, Volume 218, Issue 1, pp. 262-277, 2012. 

[15] E. Hairer, S.P. Norsett, G. Wanner, Solving Ordinary Differential Egua- 
tions 1. Nonstiff Problems. 2nd edition. Springer series in computational 
mathematics. Springer Verlag, 1993. 

[16] C.A. Kluever, Using Edelbaum’s Method to Compute Low-Thrust Trans¬ 
fers with Earth-Shadow Eclipses, Journal of Guidance, Control, and Dy¬ 
namics, Vol. 34, No. 1, January-February 2011. 

[17] H. Kuninaka, K. Nishiyama, I. Funakai, Tetsuya, Y. Shimizu, J 
Kawaguchi, Asteroid Rendezvous of HAYABUSA Explorer Using Mi¬ 
crowave Discharge Ion Engines IEPC-2005-10. Presented at the 29th 
International Electric Propulsion Conference, Princeton University, Oc¬ 
tober 31 - November 4, 2005. 

[18] W.S. Koon, M.W. Lo, J.E. Marsden, S.D. Ross, Dynamical Systems, the 
Three-Body Problem and Space Mission Design, Springer-Verlag New 
York, 2011. 

[19] G. Mingotti, F. Toputto, F. Bernelli-Zazzera, Optimal Low-Thrust In¬ 
variant Manifold Trajectories via Attainable Sets, Journal of Guidance 
Control and Dynamics, Vol. 34, No.6, 2011. 


18 



[20] G. Mingotti, P. Gurfil, Mixed Low-Thrust Invariant-Manifold Transfers 
to Distant Prograde Orbits Around Mars. Journal of Guidance, Gontrol, 
and Dynamics, Vol. 33, No. 6, November-December 2010. 

[21] G. Mingotti, F. Topputo, F. Bernelli-Zazzera, A Method to Design Sun- 
Perturbed Earth-to-Moon Low-Thrust Transfers with Ballistie Capture., 
XIX Gongresso Nazionale AIDAA, Forle, Italia, 17-21 September 2007. 

[22] G. Russell, V. Angelopoulos, The ARTEMIS Mission (Google eBook), 
Springer Science & Business Media, Nov 18, 2013 - Science - 112 pages. 

[23] T.H. Sweetser, S.B. Broschart, V. Angelopoulos, G.J. Whiffen, D.G. 
Folta, M-K. Ghung, S.J. Hatch, M.A. Woodard, ARTEMIS Mission De¬ 
sign. Space Science Reviews, Volume 165, Issue 1-4, pp. 27-57, December 
2011 . 

[24] M. Ozimek, K. Howell, Low-Thrust Transfers in the Earth-Moon System 
Ineluding Applieations to Libration Point Orbits, Journal of Guidance, 
Gontrol, and Dynamics, Vol. 33, No. 2, March-April 2010. 

[25] G. Picot, Shooting and numerieal eontinuation method for eomputing 
time-minimal and energy-minimal trajeetories in the Earth-Moon system 
using low propulsion. Discrete Gont. Dyn. Syst. Ser. B, 17, pp. 245-269, 
2012 . 

[26] L. S. Pontryagin, V. G. Boltyanskii, R. V. Gamkrelidze and E. F. 
Mishchenko, The Mathematieal Theory of Optimal Proeesses, John Wi¬ 
ley & Sons, New York, 1962. 

[27] M. Vaquero, K. Howell, Leveraging Resonant-Orbit Manifolds to Design 
Transfers Between Libration-Point Orbits, Journal of Guidance, Gon¬ 
trol, and Dynamics, Vol. 37, No. 4 (2014), pp. 1143-1157. 

[28] A. Waechter, L.T. Biegler, On the Implementation of an Interior-Point 
Pilter-Line Seareh Algorithm for Large-Seale Nonlinear Programming, 
Research Report RG 23149, IBM T.J. Watson Research Genter, York- 
town. New York, 2006. 

[29] K. Zhu, F. Jiang, J. Li, H. Baoyin, Trajectory Optimization of Multi- 
Asteroids Exploration with Low Thrust, Transactions of the Japan So¬ 
ciety for Aeronautical and Space Sciences, Vol. 54, No. 175, pp.47-54, 
2009. 


19 



